Adaptive control of CO2 bending vibration: deciphering field-system dynamics 
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We combined adaptive closed-loop optimization, phase-shaping with a restricted search space and 
imaging to control dynamics and decipher the optimal pulse. The approach was applied to controlling 
the amplitude of CO2 bending vibration during strong-field Coulomb explosion. The search space 
was constrained by expressing the spectral phase as a Taylor series, which generated pulses with 
characteristics commensurate with the natural physical features of this problem. Optimal pulses 
were obtained that enhanced bending by up to 56% relative to what is observed with comparably 
intense, transform limited pulses. We show that (f) this judicious choice of a reduced parameter 
set made unwrapping the dynamics more transparent and (2) the enhancement is consistent with 
field-induced structural changes to a bent excited state of COj^ , which theoretical simulations have 
identified as the state from which the explosion originates. 

PACS numbers: 33.15.Dj, 33.80.Eh, 42.50.Hz, 82.53.Kp 
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The confluence of intense laser fields of ultra-wide spec- 
tral bandwidth, pulse shaping [H and adaptive feedback 
[H provides a unique opportunity to explore the manipu- 
lation and engineering of many-particle dynamics at the 
quantum level. From two-particle collective motion to 
the complexities of biological interactions, the potential 
to guide an arbitrary system with tailored pulses in a 
predetermined way has opened new vistas for control via 
light-matter interaction. While much progress has been 
made towards developing schemes for realizing control - 
achieving a specific goal [1, 0] - an ability to relate the 
complex field patterns of the control pulse to a sequence 
of steps along and/or between potential surfaces describ- 
ing specific states of the system remains illusive, in gen- 
eral. Two of the primary impediments inhibiting advance 
are the vast number of (1) degrees of freedom available to 
the system under control and (2) phase and/or amplitude 
parameters influencing the pulse that must be set. To ap- 
preciate the magnitude of the latter, for phase only shap- 
ing with just 128 spectral divisions each with say 700 pos- 
sible steps between and 27r, there are ~ 1.5 x 10"^^^ per- 
mutations. Suggestions for reducing the set have ranged 
from limiting the choices to binary values (0 or tt) to 
imposing a specific functional form (see [1, 01 j ^^r ex- 
ample). In addition to reducing the search space, the 
"right" choice - one natural to the physics of the prob- 
lem - can be quite powerful, facilitating efficient genetic 
(GA) or evolutionary algorithm searches and, more im- 
portantly, enabling solutions to be deciphered physically. 
In this paper we demonstrate the power of a reduced 
set pulse shaping in a GA-mediated closed-loop control 
experiment designed to enhance bending of CO2 during 
strong-field Coulomb explosion at 800 nm. Specifically, 
we enhanced bending for the symmetric 6-electron chan- 
nel, C0^+ ^ 02+ + + 02+ , in which p\\ = for 
C2+, p|| is equal and opposite for 0^+ and p± for 0^+ is 
equal and opposite to the sum of p± for 0^+, where P\\^± 



is the momentum relative to the E field of the light. 

We chose CO2, which is linear in its ground state 
(A^f, — 0, the deviation away from 180°), as our tar- 
get because it is a nontrivial molecular system exhibiting 
modes common in larger systems ~ stretching and bend- 
ing - and it is well known to bend significantly (A^f, > 0) 
during strong- field induced Coulomb explosion d, [l^] . 
The large bending has also be the subject of theoretical 
studies, which s ugg est a physical explanation for the dis- 
tortion Ul, [ij, ll3| . The CO2 response to stron g fie lds 
allows us to rely on Coulomb explosion imaging •14 to 
probe the connection between the optimal field solutions 
and the dynamics, where the strong field pulse is the pro- 
tagonist, staring the control agent and the detector. Our 
goal was to verify the theoretical model experimentally 
by exploiting it to enhance the distortion. 

In our experiment, 50 fs transform limited (TL) pulses 
from a Tiisapphire laser system were shaped by a 128 
element phase only liquid crystal spatial light modula- 
tor (SLM). Pulses were focused with a spherical mirror 
(f.l. 75 mm) to an ~ 8 /im waist in a chamber contain- 
ing ^ 5 X 10~^ Torr of CO2. The ions generated by 
the Coulomb explosion were detected with our 47r image 
spectrometer [Tsl, [l6| . MicroChannel plates backed by a 
phosphor generate visible images of 2-dimensional pro- 
jections of 3-dimensional momentum distributions of the 
ions. Images were collected with an analog camera capa- 
ble of streaming frames to disk at 15 Hz or a fast-frame 
digital camera with a frame-storage rate of 730 Hz (see 
in Fig. [T|) . The laser rep rate was set to match the digital 
camera rate to ensure one shot per frame. Typically, the 
digital (analog) composite images contained ~ 2, 500, 000 
15, 000) laser shots. 

Searches were performed on analog images using the 
length of the 0^+ lobe as the fitness parameter. This 
was possible since 0^+ is constrained to motion _L to 
so larger A9b produces longer lobes. The search space 
was reduced by restricting the phase mask such that the 
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FIG. 1: (Color online.) Images of the symmetric 6-electron 
Coulomb explosion channel of COj"'' induced by a shaped 
pulse (top row, Fig. [21: (a) digital composite with 2,500,000 
shots; (b) selective average [1^, all frames in (a) with C^'^ 
landing inside the filled square; (c) triple-coincidence [20l |. 
all frames in (b) with ions landing simultaneously in pre- 
selected areas obeying momentum conservation (filled and 
hollow squares) and (d) analog composite of 15,000 shots. In 
all panels, E is horizontal and the center is the center of mass 
of the explosion. Panel (b) defines the far-field angle between 
the two O^""" momenta, 6cm, 147° in this case. 



spectral phase was expressed as a 5* -order Taylor ex- 
pansion p,0|, (p{uj) = I]^=o V«('^o)('^ - wo)"/n!, where 
ipni^o) = 5"</3(a;)/9w"|cjo. We ignored the first two 
terms of the series because they determine the "abso- 
lute" phase and group delay, neither of which is impor- 
tant here. Thus, the GA only varied four parameters, 
creating a phase mask by adjusting the pixels of the 
SLM collectively. It is clear from Fig. [2] that solutions 
are pulse trains that can be described in terms of even 
and odd spectral phase orders, giving us important in- 
sight into the pulse and, the dynamics it induces, as we 
discuss below. Characterization of the solution is critical 
for this study. A transie nt-g rating FROG (Frequency- 
Resolved Optical Gating 17]) was our primary tool for 
this purpose; Fig. [5] shows the FROG traces for the GA 
solutions as well. Consistency checks and calibration of 
the FROG phase retrieval were verified with a SPIDER 
(Spectral Phase Interferometry for Direct Electric-field 
Reconstruction [13]) and a SEA-TADPOLE (Spatial En- 
coded Arrangement for Temporal Analysis by Dispersing 
a Pair of Light E-fields |JJi] ) arrangements for several dif- 
ferent pulse shapes. 

Once an optimal pulse was found, statistical correla- 
tion techniques - image labeling \1M. selective averaging 
[lo| . and coincidence imaging |20|| - were run on digi- 
tal images to identify collision partners and to measure 
relative yields vs. A6b- Results are displayed in Fig. 
131 To determine A9b, we measured the far-field center 
of mass angle, 9cm, defined in Fig. [Tla, from the loca- 
tions of the correlated O^^ obeying momentum conser- 
vation (the hollow squares in Fig. [T|d). The bond angle, 
Of, = 180° — A9b, at the time of explosion was determined 
from 6cAi numerically from the equations of motion as- 
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FIG. 2: (Color online.) Reconstructed distributions from 
FROG traces (insets) for the two GA solutions (left and mid- 
dle panels) and their phase reversals (right panel) correspond- 
ing to the filled and hollow diamonds (upper row) and trian- 
gles (lower row) in Fig. |3l the solid black curves are the spec- 
tra and the temporal intensities (left axes) and the dashed 
red curves are the spectral and temporal phases (right axes). 
The reversed cases have the same spectra. The vertical axes 
of the FROG traces (insets) are the wavelength ranging from 
768 nm to 837 nm and horizontal the delay time from -1666 fs 
to 1666 fs. The parameters for the labeled peaks are given in 
Table H] Coefficients for the n — 2, ..,5 spectral phase terms 
are respectively: 800 fs^ -1.98 x 10^ fs^ -3.47 x lO** fs* and 
-5.30 X 10* fs^ (upper left panel) and -2954 fs^ -5.86 x 10^ 
fs^ 1.74 X 10^ fs'' and 1.11 X 10^ fs'^ (lower left panel). For 
the reverse cases, the sign is reversed in all terms. 

suming a pure Coulomb explosion. We point out that 
SbiScAi) depends weakly on the value of Rex [3, the 
explosion bond length, and the charge, qe, deviates by 
less than 1% for R^x in the 2.3 to 4.1 atomic units (a.u.) 
range and q between 1 to 2. Thus, any residual bond- 
ing and/or variation in R^x will have little affect on the 
angles we report. Distributions similar to that shown in 
the inset of Fig. [3] for different angles were obtained by 
selecting different C^'*' ions along the central lobe. 

Our investigation was composed of three experiments 
to probe the explosion response to the field. To isolate 
the shaped-pulse effects from changes induced by merely 
varying the intensity or the duration of the pulse, we 
performed two experiments with TL pulses. First, we 
measured the Ob distribution vs. / for 50 fs pulses, the 
solid curve in Fig. [3l The half width of the distribu- 
tion, Ai/2 db, decreases monotonically with increasing J, 
- 14° when / > 1.8 x lO^^ W/cm^, increasing to ~ 24° 
at ^ 9 X lO^'' W/cm^, the lowest intensity at which we 
could analyze images due to lack of signal strength, a 71% 
enhancement. Second, we measured the distribution for 
a 100 fs TL pulse at an intensity commensurate with the 
shaped-pulse intensity (^ 7 x 10^"^ W/cm^), the half filled 
square in Fig. [31 While the signal yield was considerably 
stronger with the longer pulse, the bending response did 
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FIG. 3: (Color online.) Bond angle distributions obtained 
from a triple-coincidence measurement and their Gaussian 
fits (upper right, for points a and 0) and vs. / for 

50 fs transform limited pulses (filled red circles), 100 fs trans- 
form limited pulse (half-filled green square) and four shaped 
pulses, the two GA solutions in the upper and lower panels 
of Fig. [2] (filled diamond and triangle respectively) and their 
reversals (hollow diamond and triangle) as discussed in the 
text. The errors in intensity reflect the fluctuations (stan- 
dard deviation) in the power and pulse-width measurements, 
while those in Ai/2^b are due to the standard deviations in 
the fit parameters. The cartoon defines the angles. 



not improve over the best 50 fs pulse response. We note 
that decreased as / increased for 100 fs pulses as 

well. Finally, we measured the distribution for two GA 
solutions, the filled triangle and diamond in Fig. [31 We 
plotted the GA solution at the intensity of the largest 
peak in the train. From the table of pulse train parame- 
ters (Table m the width of the largest peak is 100 (128) 
fs for the upper (lower) solution in Fig. [21 Clearly, the 
pulse train induces considerably more bending than does 
a TL pulse, providing an additional enhancement of 46 - 
56% (compare the 100 fs TL pulse response with that of 
the two GA solutions). 

We begin our discussion of the results by looking at 
the system response to TL fields. The sensitivity to / 
depends on several factors. First, it is well known that 
enhanced ionization 



21|, 122, [231, mediated by over-the- 



barrier ionization, is the principal ionization mechanism 
responsible for Coulomb explosion for pulses longer than 
30 fs 2J, [25|. The critical bond length, where the 
barrier is lowest, is system dependent, going as ^ 2.5//p 
a.u. for linear tiatomics j26| where Ip is the atomic ioniza- 
tion potential. For CO2, Rg is in the 3.5 to 5 a.u. range; 
we measured Rq — 4 a.u. [l^. Second, theoretical sim- 
ulations suggest the first two electrons are removed from 
CO2 before the explosion [HI Hp . which is consistent 
with our earlier measurements [10| . Third, the same sim- 
ulation shows COj"^ is promoted to a bent excited state 
prior to Coulomb explosion. The theory focused on the 
first two adiabatic states for 00^"*" (m = 0, 1,2), llVm-i- 
and \2)rn+, and found that |2)2+ was bent (Fig. 7 in |12| ) 



TABLE I: Pulse train parameters for the GA solutions in 
Fig.[5]with the (top/bottom) values corresponding to the (up- 
per/lower) GA solutions and the column headings correlate 
with the peaks in Fig. [51 n (fs) is the pulse duration; ddij 
(fs) is the period since the previous peaks; Ei is the energy of 
the i*'' peak relative to the total energy in the shaped pulse; 
and li ( xlO" W/cm^) is the intensity of the i*'' peak. The 
corresponding values for the largest peak of the reversal for 
each GA solution are; 100/104 for n, 0.62/0.59 for Si and 
7.2/6.8 for h. 
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n 100/128 


47/52 


46/48 


44/N.A. 


Mij 120/127 


87/87 


87/N.A. 


N.A./N.A. 


£i 0.63/0.57 


0.17/0.19 


0.08/0.06 


0.06/N.A. 


h 7.5/5.6 


4.2/4/3 


2.1/1.5 


1.6/N.A. 



while the other five were not. The resonant frequency 
for the G0\'^ transition, |1)2+ ^ 12)2+ , was calculated 
to close to 800 nm near Rc uM- When a system makes 
a transition from a linear state to a bent state it will be 
end up in a high vibrational level and vibrate with larger 
amplitude. Fourth, the field-dressed states are fiatter in 
the R coordinate than the field-free states (see Figs. 3, 
6 and 9 in [ll[) and the dressed |2)2+ state is less bent, 
A^f, ~J.5°, than it is field- free, A^b ~ 60° (compare Fig. 



7 in [1J| with Fig. 17 in [13|). Flatter potentials allow 
the system to stretch further making it easier to reach 
Rc- Consequently, theory suggests that Coulomb explo- 
sions originating from |2)2+ will exhibit larger bending 
amplitude than those originating from 1 1)2-1- • We expect 
less bending from |2)2+ when dressed with a stronger 
field because the |1)2+ 1 2) 2-1- transition could promote 
population to lower vibrational states in |2)2+. This par- 
tially explains our observations. We point out, however, 
that with the observed increase in the bending for de- 
creasing / (Fig. [3]), is a concomitant, monotonic decrease 
in kinetic energy release during the explosion, indicating 
Rex oc 1//. This well known saturation is a result of the 
threshold / for over-the-barrier ionization increasing as 
R decreases. For strong enough fields, over-the-barrier 
ionization and the explosion occur at R^x < Rc- Thus, 
near Rc Coulomb explosion is induced at lower /, which 
again leads a larger bending for smaller /. 

Turning to the GA solutions (Fig. [2]), we notice that 
while the two pulse trains are not identical (see Table 
[J for a summary of their characteristics), they do share 
two important features: (1) the last peak is larger and 
wider than the earlier peaks - a pulse width > 100 fs 
compared with ~ 50 fs for early pulses - and (2) the 
spectral phase is dominated by a negative odd (largely 
third) order (see the caption of Fig.[2lfor the ipn values). 
A pulse train preceding a larger peak is characteristic of 
a negative odd order spectral phase. In fact, all solutions 
shared these two features. All solutions showed bending 
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enhancement over TL pulses, albeit, with some inducing 
significantly more bending than others. One advantage 
of this reduced parameter set is that we can change the 
sign of or turn off specific phase terms. To verify that 
the negative odd order spectral phase was in fact impor- 
tant, we reversed the sign of the ipn terms to produce a 
pulse train with the first peak being largest (right pan- 
els in Fig. [21). The result was the enhancement virtually 
disappeared as shown by the hollow diamond and trian- 
gle in Fig. [3l The most likely reason the train of pulses 
enhances bending vibration is a combination of the fact 
that the |2)2+ state is populated by the early peaks more 
completely than with a single pulse and the field oscillates 
off and back on. As we pointed out earlier, atomic ions 
were essentially immeasurable for / < 7 x lO^"* W/cm^. 
Even though Coulomb explosions were insignificant, this 
intensity was sufficient to produce both COj and COj^ 
ions. So, each of the early peaks help to populate the 
12)2+ state. Simulations show that subsequent pea ks do 
not transfer population back to the |1)2+ state [27[. Be- 
tween pulses, when the field is off, the state relaxes to its 
more bent field-free condition causing CO2+ to vibrate 
more. It is also interesting to note that the period be- 
tween pulses is near the bending vibrational frequency 
of the ground state. Thus, the field could also provide a 
periodic kick, further enhancing the vibration. It is im- 
portant to see that the GA enhancement does not come 
from a lowering of /. The 100 fs TL and the two re- 
versal solutions all show significantly less bending than 
the GA solutions but have / about the same as that of 
the diamond solution. The kinetic energy release for the 
filled triangle and filled diamond are identical. The fact 
that the kinetic energy release for the GA solutions was 
identical suggest that any enhancement due to intensity 
and bond length at the time of explosion is probably 
minimal. We then conclude that the pulse train with a 
negative third (odd) order chirp plays a primary roll in 
enhancing the bending amplitude. 

The results presented here coupled with our earlier 
studies [lo'l are consistent with the theoretical model of 
the CO2 explosion originating from a high vibrational 
level of a bent state in C02^. For TL pulses, the bend- 
ing amplitude decreases with increasing /, in response to 
a straightening of the bent state for larger J. A shaped 
pulse composed of a pulse train provides additional bend- 
ing enhancement. The dressed-bent state is allowed to re- 
lax (become more bent) between peaks within the train 
amplifying the vibration, with the early peaks populat- 
ing the state and final peak inducing the explosion. A 
judicious choice for the search space proved to be power- 
ful in deciphering the optimal pulse. Vibration was the 
natural physical feature for this problem and restricting 
the phase mask to producing pulses forced the optimal 
pulse to reflect the primary dynamics. 
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